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Abstract: The understanding of the particle spectra resulting from acceleration in relativistic shocks as 
they occur in extragalactic sources, is essential for the interpretation of the cosmic ray spectrum above the 
ankle (E p > 3 ■ 10 8 eV). It is believed that extragalactic sources like Active Galactic Nuclei and Gamma 
Ray Bursts can produce particle spectra up to E v ~ 10 21 eV. In this contribution, subluminal shocks are 
investigated with respect to different shock boost factors V and the inclination angle between the shock 
normal and the magnetic field tp. A correlation between the boost factor and the spectral behavior of the 
emitted particles is found. The results are compared to Active Galactic Nuclei and Gamma Ray Burst 
diffuse cosmic ray contribution and the observed cosmic ray spectrum at the highest energies. 
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Introduction 

The observations of electron synchrotron radiation 
in the radio regime and the X-ray variations from 
astrophysical environments such as Super Novae 
(SN), Active Galactic Nuclei (AGN) and Gamma 
Ray Bursts (GRBs), indicate that there is a mecha- 
nism with which primary particles (protons or elec- 
trons) gain large amounts of energy via accelera- 
tion. Connected to this implied mechanism, it is 
accepted that shocks are the source of the accel- 
eration of the observed cosmic rays. The work 
of the late 70s by a number of authors such as 
[12, 3, 4, 1, 5] established the basic mechanism 
of particle diffusive acceleration in non-relativistic 
shocks. Later on studies on relativistic shocks ap- 
peared (e.g. [17, 11, 19, 7]) and questioned the 
nature and the limits of the relativistic diffusive ac- 
celeration versus the non-relativistic shock mecha- 
nism. In this work, we present results on relativis- 
tic particle acceleration which are extended to rela- 
tivistic source models (i.e. AGN, GRBs). These in- 
volve subluminal shocks with plasma shock veloc- 
ity ranging from Lorentz boost factors, T = 1 — > 
1000, using pitch angle and magnetic scattering. 
In the first section, the Monte Carlo simulations 
are described, while in the second section the re- 



sults are presented for different boost factors. In 
the third section, the Monte Carlo spectra are used 
to calculate the possible contribution of AGN and 
GRBs to the observed spectrum of charged cosmic 
rays. The final section summarizes the results. 

Numerical simulations 

The Monte Carlo codes simulate relativistic 
oblique shocks with application to relativistic en- 
vironments of AGN and GRBs. Initially the par- 
ticles are considered relativistic, isotropically in- 
jected upstream of the shock and they are let to 
move towards the shock and along the way they 
collide with the scattering centers. Consequently, 
as they keep scattering between the upstream and 
downstream regions of the shock (its width is much 
smaller than the particle's gyroradius), they gain a 
considerable amount of energy. During our simula- 
tions, Lorentz relativistic transformations are per- 
formed from the local plasma frames to the normal 
shock frame to check for shock crossings. 
Particles leave the system if they escape far down- 
stream at a spatial boundary r&, or if they reach a 
defined maximum energy E max , even though other 
physics describing particle escape or energy loss 
would probably need to be taken into account in re- 
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alistic situations (forthcoming paper). The down- 
stream spatial boundary required can be estimated 
from the solution of the convection-diffusion equa- 
tion in a non-relativistic, large-angle scattering ap- 
proximation in the downstream plasma. 
The mean free path is calculated in the respective 
fluid frames by the formula A — X ■ pi j2 , assum- 
ing a momentum dependence to this mean free path 
for scattering along the field, related to the spatial 
diffusion coefficient, k. In the shock normal, or in 
the x direction, the diffusion coefficient is given by 
«ll = Av/3 where k = n^cos 2 ip. At the scattering 
centers the energy (momentum) of the particle is 
kept constant and only the direction of the velocity 
vector v is randomized, using a computational ran- 
dom number generator. Runs were performed with 
different spatial boundaries to investigate the effect 
of the size of the acceleration region on the spec- 
trum, as well as to find a region where the spectrum 
is size independent. In the small-angle scattering 
case, the inherent anisotropy due to the high down- 
stream sweeping effect, may greatly modify this 
analytical estimate. A splitting technique is used 
similar to the one used in the Monte Carlo simu- 
lations of [15], [16], so that when an energy level 
is reached such that only few accelerated cosmic 
rays remain, each particle is replaced by a num- 
ber of N particles of statistical weight 1/N, so as 
to keep roughly the same number of cosmic rays 
being followed. Workers on the field [20] have 
demonstrated that small angle scattering (pitch an- 
gle diffusion), with SO < 1/T, applies with 9 mea- 
sured in the upstream fluid frame for scattering in 
a uniform field or a randomly oriented set of uni- 
form field cells. This arises because particles at- 
tempting to penetrate upstream from the shock are 
swept back into the shock before they can move 
far in pitch. For this case, we allow particles' pitch 
angle diffusion to angles chosen at random, up to 
an angle ~ 10/T, where T is the upstream gamma, 
measured in the normal shock frame. 
For the oblique shock cases that we apply here, 
provided the field directions encountered are rea- 
sonably isotropic in the shock frame, we know that 
tan tpi = r^ 1 tan -0s ~ T^ 1 ~ ipi where '1' and 
V refer to the upstream and normal shock frames 
respectively. Given the current interest in follow- 
ing shock acceleration in the relativistic flow in test 
particle regime, we outline a typical approach to 



such computations. Using the guiding center ap- 
proximation, a test particle moving a distance, d, 
along a field line at tp to the shock normal, in the 
plasma frame has a probability of collision within 
d given by P(d) = 1 - cxp(-d/A) = R. Where 
the random number R is < R < 1. Weight- 
ing the probability by the current in the field direc- 
tion, /z, then d = -A/ilni?. Between the 'i th' 
and 'i+1 th' scattering, a shock crossing condition 
is checked, by finding the position and time in the 
normal shock frame, as 

Vl/Il 

At: +l = ^(1^1^! + ^) 

where ' 1 ' refers to the quantities that are calculated 
in the upstream fluid frame. For further details 
on relativistic transformation calculations see [14]. 
We note that for the description of the physical 
quantities necessary throughout the simulations we 
use the fluid rest frames, the normal shock frame 
and the de Hoffman-Teller frame where E=0. Fur- 
thermore, as the theory of the oblique shocks pos- 
tulates, [9], from the conservation of the first adia- 
batic invariant, one can find the new pitch angle in 
the downstream frame and similar transformations 
allow the particle scattering to be followed in this 
frame. 

Results 
Energy spectra 

As discussed above, Monte Carlo calculations 
have been performed for different parameter set- 
tings. Three different shock inclinations have 
been simulated, i. e. tp = 23°, 33° and 43°, 
where tp is the angle between the shock nor- 
mal and the magnetic field. For each of the 
angles, nine different boost factors have been 
used, T = 10, 20, 30, 100, 300, 500, 700, 900 
and 1000. Results will be presented in detail 
in [14]. Furthermore, each single spectrum is fit 
by a single power-law approximation, dN/dE oc 
E~ ap , in order to get an estimate of the spectral 
index a p . At the lowest energies, particle injec- 
tion features are excluded. Since AGN and GRBs 
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Figure 1 : Correlation between the shock boost fac- 
tor r and the spectral index, for shock inclination 
of ijj = 43° . The correlation appears to be linear. 



are likely to accelerate particles up to ~ 10 21 eV, 
this value is used as the maximum energy. This 
way, each shock boost factor T is associated with a 
spectral index a p (T). A linear correlation between 
a p and T is found as it is displayed in Fig. 1. The 
figure shows the results for rp = 43°, the results 
for ip = 23° and tp = 33° are comparable. A lin- 
ear fit is applied with a numerical result which is 
comparable for all three angles, 

a p (tp = 23°) = 2.1 - 0.7- 1CT 3 T 
a p (tp = 33°) = 2.0-0.9- 1CT 3 T 
a p (ip = 43°) = 2.2-0.9-10~ 3 -r. 

This implies three main things for relativistic, sub- 
luminal shocks: 

1 . The result is basically independent of the in- 
clination angle between shock normal and 
magnetic field, as long as the shocks are 
oblique. 

2. Mildly-relativistic shocks have spectral in- 
dices of around a p ~ 2.0 — 2.4. 

3. The more relativistic the shock, the flat- 
ter the spectrum. Highly-relativistic shocks 
with T ~ 1000 can produce spectra as flat as 
a v ~ 1.5. 



The above findings seem to support similar simula- 
tion results and observational evidences regarding 
irregular and fiat spectra from many GRBs (further 
details in [14]). 

Contribution to the charged Cosmic Ray 
spectrum 

In this section, the spectra of different boost fac- 
tors are associated to different source classes in 
order to calculate the possible contribution to the 
diffuse spectrum of charged cosmic rays. Specif- 
ically, it is assumed that the cosmic ray spectrum 
above the ankle, E > 10 185 eV, is produced by 
protons accelerated in AGN and GRBs sources. 
For AGN, a fix boost factor of T = 10 is used. 
In the case of GRBs, it is assumed that 10% of all 
GRBs have T = 100, 80% have T = 300 and the 
remaining 10% have T = 1000. It is further as- 
sumed that the observed spectrum of charged cos- 
mic rays is produced by a linear combination of 
AGN, dN AGN /dE, and GRBs, dN GRB /dE, 

dN C R dN AGN dN G RB 

= a • h 1 — a • , 

dE dE y ' dE ' 

with < a < 1. Source evolution has been taken 
into account by using the source density given by 
the observation of X-ray AGN [8]. This source 
evolution is also used in the case of GRBs, assum- 
ing that GRBs also follow the Star Formation Rate. 
The energy density of the observed comic ray spec- 
trum je is then used to normalize the calculated 
Monte Carlo spectra 

JE ■= / 7F EdE 
w lO^GeVcm^s^sr- 1 . 

The final result is shown in Fig. 2. The upper line 
represents a pure AGN spectrum, a = 1, the mid- 
dle line represents a mixture of 50% AGN and 50% 
GRBs, a = 0.5, and the lower line is a pure GRB 
spectrum, a = 0. Since the cosmic ray spectrum 
above E = 10 18 5 eV needs to be explained by ex- 
tragalactic sources, the best fit results from a pure 
AGN spectrum. With a significant contribution, 
the spectrum becomes too flat. 
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Figure 2: The diffuse energy spectrum of relativis- 
tic, subluminal sources, compared to the total dif- 
fuse spectrum. 



Conclusions 

In this work we presented Monte Carlo simulation 
studies on the acceleration of test particles in rel- 
ativistic, subluminal shock environments. Source 
candidates were applied, such as Active Galac- 
tic Nuclei with mild relativistic shocks of boost 
factors of T w 10 and Gamma Ray Bursts with 
highly-relativistic shocks, 100 < V < 1000. 
The resulting particle spectra from those shocks 
have been investigated in the context of vary- 
ing shock boost factor V versus the shock obliq- 
uity, i. e. the inclination angle between the shock 
normal and the magnetic field, ip. For three 
angles, ijj = 23°, 33°, 43°, shocks were sim- 
ulated with a series of different boost factors, 
T = 10, 20, 30, 100, 300, 500, 700, 900, 1000. 
While the result remained constant for the dif- 
ferent inclination angles, a linear behavior be- 
tween the spectral index a p and the boost factor 
r was found, leading to spectra of a p ~ 2.2 for 
mildly-relativistic shocks (T ~ 10), but producing 
much harder spectra for highly-relativistic shocks. 
Specifically, for boost factors of T ~ 300 - 1000, 
the spectra can become as flat as a p ~ 1.5. 
Concerning the observed cosmic ray spectrum at 
the highest energies, our study shows that one can- 
not explain those high energies using solely a con- 



tribution from GRBs. On the other hand, a pure 
AGN spectrum fits the spectrum well. 
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